What is claimed is: 

1. A method of estimating a set of floating ambiguities associated with a set of 

phase measurements of a plurality n of satellite carrier signals made by a first navigation 
receiver (B) and a second navigation receiver (R) separated by a distance, wherein a 
baseline vector (x°,y°,z°) relates the position of the second receiver to the first receiver, 
5 each satellite carrier signal being transmitted by a satellite and having a wavelength, 

wherein each receiver has a time clock for referencing its measurements and wherein any 
difference between the time clocks may be represented by an offset, said method 
receiving, for a plurality of two or more time moments j 9 the following inputs: 

a vector representative of a plurality of pseudo-ranges measured by the first 
10 navigation receiver (B) and corresponding to the plurality of satellite carrier signals, 

a vector yf* representative of a plurality of pseudo-ranges measured by the second 
navigation receiver (R) and corresponding to the plurality of satellite carrier signals, 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the first navigation receiver (B), 
15 a vector Dj representative of a plurality of estimated distances between the 

satellites and the second navigation receiver (R), 

a vector <Pj representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the first navigation receiver (B), 

a vector q>j representative of a plurality of flill phase measurements of the 
20 satellite carrier signals measured by the second navigation receiver (R), 

a geometric Jacobian matrix Hj* whose matrix elements are representative of the 
changes in the distances between the satellites and one of the receivers that would be 
caused by changes in that receiver's position and time clock offset, said method 
comprising the steps of: 

25 (a) generating, for each time moment j 9 a vector A# of a plurality of range 

residuals of pseudo-range measurements made by the first and second navigation 
receivers in the form of: A# = (yf -jf)- (Dj R - D* ); 
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(b) generating, for each time moment j 9 a vector A <Pj of a plurality of phase 
residuals of full phase measurements made by the first and second navigation receivers in 
the form of: 

= (fc R - qf) - A"' ■(/>/ - Df ), 

where A is a diagonal matrix comprising the inverse wavelengths of the satellites; 

(c) generating, for time moment j — 1, an LU- factorization of a matrix Mj or a 

matrix inverse of matrix Mi, the matrix Mi being a function of at least A and Hp\ 

(d) generating, for time moment j = 1, a vector Ni as a function of at least A?'/, 

A g>j 9 and the LU- factorization of matrix Mi or the matrix inverse of matrix Mi; 

(e) generating, for an additional time moment j 5*1, an LU-factorization of a 

. -l 

matrix My or a matrix inverse of matrix My, the matrix My being a function of at least A , 
Hp and an instance of matrix M generated for a different time moment; and 

(f) generating, for an additional time moment j 5*1, a vector Ny as a function of at 
least A75, A q>j , and the LU-factorization or matrix My or the matrix inverse of matrix My, 
the vector Ny having estimates of the floating ambiguities. 

2. The method of Claim 1 wherein step (c) comprises generating an LU- 



• the matrix G has 2n rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprises annxn zero matrix and the other sub-matrix 
comprising annxn identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

• the matrix Pi has 2n rows, 2n columns, and a form which comprises a matrix 

equivalent to ^Rf 1 - Rf 1 Qi(q? Rf 1 Qi +^S 1 jT 1 Qj r Rf 1 j where the matrix Ri 

is a weighting matrix, where the matrix R1" 1 comprises an inverse of matrix Ri, where 
the matrix Qi has an upper sub-matrix and a lower sub-matrix, one of the sub- 
matrices of Qi comprising matrix Hj y and the other of the sub-matrices of Qi 

comprising the matrix product A Hp 9 and wherein the matrix Qi T comprises the 
transpose of matrix Qi, and where the quantity qSk is a zero matrix when the distance 



factorization for a matrix comprising a form equivalent to (G PjG I, where: 
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between the first and second navigation receivers is unconstrained and where q may 
be a non-zero weighting parameter and may be a non-zero matrix when the 
distance between the first and second navigation receivers is constrained. 

3. The method of Claim 2 wherein step (d) comprises the step of generating 

vector Ni to comprise a vector having a form equivalent to M j" 1 (g T P^ fij + qgi ), where 
the matrix Mi" 1 comprises an inverse of matrix of matrix Mi, and where the vector [i\ 
comprises the vector [A/i, Aq>\ ] T , and where the quantity qgk is a zero vector when the 
distance between the first and second navigation receivers is unconstrained and where q 
may be non-zero and gk may be a non-zero vector when the distance between the first and 
second navigation receivers is constrained. 

4. The method of Claim 1 wherein step (e) comprises generating an LU- 

T 

factorization for a matrix comprising a form equivalent to My = My_j + G P/G 9 
where: 

• Mj-j comprises the matrix Mi of step (c) when j = 2 and comprises the matrix M y of 
step (e) for the j-l time moment when j > 2, 

• the matrix G has 2n rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprising annx/i zero matrix and the other sub-matrix 
comprising annxn identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

• the matrix P, has In rows, In columns, and a form which comprise a matrix 

equivalent to ^R} 1 - Ry 1 Q 7 (qJ RJ 1 Qj + ?S 7 )^ 1 qJ RJ*j where the matrix 

Ry is a weighting matrix, where the matrix R/~ l comprises an inverse of matrix Ry, 
where the matrix Q y has an upper sub-matrix and a lower sub-matrix, one of the sub- 
matrices of Qj comprising matrix Hj* and the other of the sub-matrices of Q, 

comprising the matrix product A H/ 9 and wherein the matrix Q 7 T comprises the 
transpose of matrix Q/, and where the quantity qS y is a zero matrix when the distance 
between the first and second navigation receivers is unconstrained and where q may 
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be a non-zero weighting parameter and S, may be a non-zero matrix when the distance 
between the first and second navigation receivers is constrained. 

20 

5. The method of Claim 4 wherein step (f) comprises the step of generating vector 
Nj to comprise a vector having a form equivalent to 

Ny_i -f My 1 [g T P/ ( li j - GNy_i )+ qgj ], where the matrix M/ 1 comprises an 

inverse of matrix of matrix M y , where the vector fij comprises the vector [A#, A <Pj] T , and 
5 where the vector Ny_/ comprises the vector Ni generated by step (d) when j = 2 and 

comprises the vector N/_; generated by step (f) for the j-l time moment when j > 2, and 
where the quantity qgj is a zero vector when the distance between the first and second 
navigation receivers is unconstrained and where q may be non-zero and gj may be a non- 
zero vector when the distance between the first and second navigation receivers is 
10 constrained. 



10 



6. The method of claim 3 wherein the distance between the first and second 
navigation receivers is constrained to a distance Lrb, wherein step (c) generates matrix Si 
in a form equivalent to: 



Si = 



V 



1- 



'1x3 



°3xl i , Lrb T 

o J + M r ' ri 



where Y\ is a vector comprising estimates of the three coordinates of the baseline vector 
for the time moment j=l and a zero as fourth component, where Ti T is the vector 
transpose of where I 3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three 
zeros, and where C>3 X i is a column vector of three zeros; and 

wherein step (d) generates vector gi for the time moment j=l in a form equivalent 

to: 



gl =G T Rr 1 Qi(Q 1 T Rf 1 Qi+qrS 1 r 1 h 1 , 



where: 



hi = 



1 L 

1 UP 
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7. The method of Claim 3 wherein weighting matrix Ri comprises an identity 
matrix multiplied by a scalar quantity. 



8. The method of claim 5 wherein the distance between the first and second 
navigation receivers is constrained to a distance Lrb, wherein step (e) generates matrix S,- 
in a form equivalent to: 

( , V 



1- 



JWJ 



{ i 3 o 3xl ~ 

^>i x3 0 , 



'A 



where r,- is a vector comprising estimates of the three coordinates of the baseline vector 
for the j-th time moment and a zero as fourth vector component, where r 7 T is the vector 
transpose of r,-, where I3 is the 3-by-3 identity matrix, where Oj X 3 is a row vector of three 
zeros, and where 0 3x i is a column vector of three zeros; and 

wherein step (f) generates vector g/ for the j-th time moment in a form equivalent 

to: 

gj =G T R-%(Q j T RJ 1 Q J + q S j r 1 h j , 

where: 



h 7 = 



1 



J\\J 



9. The method of Claim 5 wherein the weighting matrix R, comprises an identity 
matrix multiplied by a scalar quantity for at least one time moment j. 

10. The method of Claim 4 wherein the generation of the LU- factorization in step 
(e) comprises the steps of: 

(g) generating an LU-factorization of matrix Mj.j in a form equivalent to 
Ly_i Ltj-\ T wherein L/_i is a low-triangular matrix and L/_i r is the transpose of L 7 _i; 

(h) generating a factorization of G T P y G in a form equivalent to T, T y T = G T P, G, 
where T/ is the transpose of T y ; 
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(i) generating an LU- factorization of matrix My in a form equivalent to L, L y r from 
a plurality n of rank-one modifications of matrix Ly_i, each rank-one modification being 
based on a respective column of matrix Ty, where n is the number of rows in matrix My . 

11. The method of Claim 10 wherein step (h) generates matrix Ty from a 
Cholesky factorization of G T Py G. 



12. The method of Claim 10 wherein weighting matrix Ry has a form equivalent 



to: 



R j = 



R) 0 
0 Rj 



, where R Y and R^ are weighting matrices; 



wherein R Y and R^ 7 are related to a common weighting matrix W and scaling 
parameters oy and as follows (r y )~ = — jW, and (r^)^ ^^-W, 

wherein step (h) of generating matrix Ty comprises the steps of: 

(Ty A*QPS 

generating a scalar b in a form equivalent to: b = — 2 2 2 9 w ^ ere ^GPS * s 

°> ~ h ^GPS <J <p 

the wavelength of the satellite signals, 

generating a matrix H in a form equivalent to H = W 1/2 H^ , 

generating a Householder matrix S H h for matrix H , and 
generating matrix Ty in a form equivalent to: 



Ty =— W 1/2 Shh 



V, A GPS J 



p 4x(«-4) 



0(„_4) x4 | I( rt -4) x («-4) 



13. The method of Claim 10 wherein weighting matrix Ry is applied to a case 
where there is a first group of satellite signals having carrier signals in a first wavelength 
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band and a second group of satellite signals having a carrier signals in a second 
wavelength band, the weighting frequency having a form equivalent to: 



R) 0 
0 Rj 



, where R Y and are weighting matrices; 



: earner 



wherein R Y and R^ are related to a common weighting matrix W, the carrier 

wavelengths of the first group of signals as represented by matrix , the < 

wavelengths of the second group of signals as represented by matrix A< 2 > , the center 
wavelength of the first band as represented by , the center wavelength of the first band 

as represented by A 2 , and scaling parameters a y and cr^ , as follows: 

1 



tor- 



w 



o 



nxn 



o 



nxn 



w 



kr - 



>q> 

o 



w 



(1) 



o 



nxn 



nxn 



w 



(2) 



15 



where 



wherein step (h) of generating matrix T) comprises the steps of: 



generating a scalar b in a form equivalent to: b — 



where is the wavelength of a first group of satellite signals and X 2 is the wavelength 
of a second group of satellite signals, 

generating a matrix H in a form equivalent to H = W 1 1 2 H^ , 
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20 



25 



generating a Householder matrix Shh for matrix H , and 



generating matrix T, in a form equivalent to: Tj = 



All | O nxn 
A21 | A22 



where sub-matrixces All, A21, and A22 are as follows: 

^\-b/Xx I 4 I 0 4x(n _ 4) 



All 



= (W< 1 ))2 ! 



»HH 



°(n-4)x4 I I (n-4) 



A21 



= (w< 2 >)* 



2 S 



HH 



l 4 I °4x(n-4) 



°(n-4)x4 I °(n-4)x(n-4) 



, and 



A22 



= (\V< 2 >)2| 



'HH 



\-b/xf -b/xt; 
i-b/xf 



4x(n-4) 



I 

°(n-4)x4 I I (n-4) 



10 



14. The method of Claim 4 wherein step (d) comprises generating a vector Bi to 
comprise a vector having a form equivalent to G T Pi J*i+qgi, where the vector fL\ 
comprises the vector [Ayi, dq>\] T , and where the quantity qgj is a zero vector when the 
distance between the first and second navigation receivers is unconstrained and where q 
may be non-zero and gj may be a non-zero vector when the distance between the first and 
second navigation receivers is constrained; and 

wherein step (f) further comprises generating, for each time moment j 5*1, a vector 
Bj to comprise a matrix having a form equivalent to B^i + G T P y fij +qg/, where the vector 
/t, comprises the vector [A75, A <pj ] T , and where the vector B/_/ is the vector Bi generated 
by step (d) when j = 2 and comprises the vector generated by step (f) for the for the j-l 
time moment when j > 2, and where the quantity qgy is a zero vector when the distance 
between the first and second navigation receivers is unconstrained and where q may be 



66 



non-zero and g, may be a non-zero vector when the distance between the first and second 
navigation receivers is constrained; and 

wherein step (f) further comprises generating vector N, to comprise a vector 
having a form equivalent to N y = [M/]~ l By, where the matrix M/ 1 comprises an inverse of 
matrix of matrix My. 



15. The method of claim 14 wherein the distance between the first and second 
navigation receivers is constrained to a distance Lrb, wherein step (c) generates matrix Si 
in a form equivalent to: 



Si = 



1- 



RB 

NU 



v 



'3x1 



^1x3 



+ ^r 1 ri 7 ' 



where Y\ is a vector comprising estimates of the three coordinates of the baseline vector 
for the time moment j=l and a zero as fourth component, where i*i T is the vector 
transpose of where I3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three 
zeros, and where 0 3x i is a column vector of three zeros; and 

wherein step (d) generates vector gj for the time moment j=l in a form equivalent 

to: 



where: 



1- 



L RB 

Ikill 



16. The method of claim 14 wherein the distance between the first and second 
navigation receivers is constrained to a distance Lrb, wherein step (e) generates matrix S, 
in a form equivalent to: 



Sj = 



1- 



L RB 



i 3 o 3xl 
Lo lx3 0 



where r, is a vector comprising estimates of the three coordinates of the baseline vector 
for the j-th time moment and a zero as fourth vector component, where r 7 T is the vector 
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transpose of r,, where I3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three 
zeros, and where O3X1 is a column vector of three zeros; and 

wherein step (f) generates vector g, for the j-th time moment in a form equivalent 

to: 

gj =g t rj 1 q 7 (q/r7 1 q 7 +< 7 s y r 1 h 7 , 

where: 

f \ 

— ft- 



J 




1 7. The method of Claim 14 wherein the weighting matrix R, comprises an 
identity matrix multiplied by a scalar quantity for at least one time moment j. 



18. A method of estimating a set of floating ambiguities associated with a set of 
phase measurements of a plurality n of satellite carrier signals made by a first navigation 
receiver (B) and a second navigation receiver (R), wherein a baseline vector (x°,y°,z°) 
relates the position of the second receiver to the first receiver, each satellite carrier signal 
being transmitted by a satellite and having a wavelength, wherein each receiver has a time 
clock for referencing its measurements and wherein any difference between the time 
clocks may be represented by an offset, said method receiving, for a plurality of two or 
more time moments j, the following inputs for each time moment j: 

a vector yf representative of a plurality of pseudo-ranges measured by the first 
navigation receiver (B) and corresponding to the plurality of satellite carrier signals, 

a vector yf~ representative of a plurality of pseudo-ranges measured by the second 
navigation receiver (R) and corresponding to the plurality of satellite carrier signals, 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the first navigation receiver (B), 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the second navigation receiver (R), 

a vector q)j representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the first navigation receiver (B), 
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R 

a vector q>j representative of a plurality of full phase measurements of the 
20 satellite carrier signals measured by the second navigation receiver (R), 

a geometric Jacobian matrix Hj whose matrix elements are representative of the 
changes in the distances between the satellites and one of the receivers that would be 
caused by changes in that receiver's position and time clock offset, said method 
comprising the steps of: 

25 (a) generating, for each time moment j 9 a vector A# of a plurality of range 

residuals of pseudo-range measurements made by the first and second navigation 
receivers in the form of: A/j= (jf - }f) - (D* -/>/), said step generating a set of 
range residuals A% , k=l 9 ..,7; 

(b) generating, for each time moment j, a vector A <Pj of a plurality of phase 

30 residuals of full phase measurements made by the first and second navigation receivers in 
the form of: A = (q>j - tpj ) — A (Dj - Dj ), where A is a diagonal matrix 
comprising the inverse wavelengths of the satellites, said step generating a set of phase 
residuals A , fc=l , ..,y; 

(c) generating an LU-factorization of a matrix M or a matrix inverse of matrix M, 
35 the matrix M being a function of at least A and //* 7 , for index k of H k y covering at 

least two of the time moments j ; 

(d) generating a vector N of estimated floating ambiguities as a function of at least 
the set of range residuals A% , the set of phase residuals A q>^ and the LU-factorization of 
matrix M or the matrix inverse of matrix M. 



19. The method of Claim 18 wherein step (c) comprises generating matrix M in a 



form equivalent to the summation 



£(G T P*G) 

k=\ 



where: 

• the matrix G has In rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprises annxn zero matrix and the other sub-matrix 
comprising annxn identity matrix, 



• the matrix G comprises the transpose matrix of matrix G, and 
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• the matrix P* has 2n rows, 2n columns, and a form which comprises a matrix 

equivalent to = R* 1 - R* 1 Q k {Qj R* 1 Q k + qS k )~ l Qj R* 1 , where the 
10 matrix R^ is a weighting matrix, where the matrix R*" 1 comprises an inverse of matrix 

R*, where the matrix has an upper sub-matrix and a lower sub-matrix, one of the 
sub-matrices of comprising matrix H k and the other of the sub-matrices of 
comprising the matrix product A //* 7 , and wherein the matrix Q* T comprises the 
transpose of matrix Q*, and where the quantity qSk is a zero matrix when the distance 
1 5 between the first and second navigation receivers is unconstrained and where q may 

be a non-zero weighting parameter and S* may be a non-zero matrix when the distance 
between the first and second navigation receivers is constrained. 



20. The method of Claim 19 wherein step (d) comprises generating matrix N in a 
form equivalent to: 



N = M _1 x 



£(G T P*^+<7g k ) 

k=\ 



where the matrix IVT 1 comprises an inverse of matrix of matrix M, where the vector n k 
comprises the vector [Ay k , Aq> k ] T , and where the quantity qgk is a zero vector when the 
distance between the first and second navigation receivers is unconstrained and where q 
may be non-zero and gk may be a non-zero vector when the distance between the first and 
second navigation receivers is constrained. 



21 . The method of claim 20 wherein the distance between the first and second 
navigation receivers is constrained to a distance Lrb, wherein step (c) generates matrix Sk 
for the k-th time moment in a form equivalent to: 



f 



1- 



13 o 3xl 

Pi x3 0 



, L RB „ T 



where r* is a vector comprising estimates of the three coordinates of the baseline vector at 
the k-th time moment, and a zero as fourth component, where r* T is the vector transpose 
of I*, where I3 is the 3-by-3 identity matrix, where Oj X 3 is a row vector of three zeros, and 
where 0 3x i is a column vector of three zeros; and 
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wherein step (d) generates vector gk for the k-th time moment in a form equivalent 

10 to: 

g* = G T R* 1 Q*(Q* T R* 1 Q* +gS k y 1 h k , 

where: 

( T ^ 
I IF* I J 

22. The method of Claim 19 wherein at least one of the weighting matrices Rk 
comprises an identity matrix multiplied by a scalar quantity. 

23, A computer program product for directing a data processor to estimate a set of 
floating ambiguities associated with a set of phase measurements of a plurality n of 
satellite carrier signals made by a first navigation receiver (B) and a second navigation 
receiver (R) separated by a distance, wherein a baseline vector (x°,y°,z°) relates the 

5 position of the second receiver to the first receiver, each satellite carrier signal being 

transmitted by a satellite and having a wavelength, wherein each receiver has a time clock 
for referencing its measurements and wherein any difference between the time clocks may 
be represented by an offset, the process receiving, for a plurality of two or more time 
moments j 9 the following inputs: 

10 a vector yf representative of a plurality of pseudo-ranges measured by the first 

navigation receiver (B) and corresponding to the plurality of satellite carrier signals, 

a vector yf representative of a plurality of pseudo-ranges measured by the second 
navigation receiver (R) and corresponding to the plurality of satellite carrier signals, 

a vector Dj representative of a plurality of estimated distances between the 
15 satellites and the first navigation receiver (B), 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the second navigation receiver (R), 

a vector representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the first navigation receiver (B), 



71 



a vector <pj representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the second navigation receiver (R), 

a geometric Jacobian matrix Hj* whose matrix elements are representative of the 
changes in the distances between the satellites and one of the receivers that would be 
caused by changes in that receiver's position and time clock offset, the computer program 
product comprising: 

a computer-readable medium; 

a first set of instructions embodied on the computer-readable medium which 
directs the data processor to generate, for each time moment j, a vector AjJ of a plurality 
of range residuals of pseudo-range measurements made by the first and second navigation 
receivers in the form of: = ( yf^ — yf) — (Dj R — Df ); 

a second set of instructions embodied on the computer-readable medium which 
directs the data processor to generate, for each time moment j 9 a vector A q)j of a plurality 
of phase residuals of full phase measurements made by the first and second navigation 

receivers in the form of: A <Pj = (<Pj K - ^ B ) — A (Pf - Dj B ), where A is a diagonal 
matrix comprising the inverse wavelengths of the satellites; 

a third set of instructions embodied on the computer-readable medium which 
directs the data processor to generate, for time moment j = 1, an LU- factorization of a 
matrix Mi or a matrix inverse of matrix Mi, the matrix Mi being a function of at least 

A~ l and Hj y ; 

a fourth set of instructions embodied on the computer-readable medium which 
directs the data processor to generate, for time moment j = 1, a vector Ni as a function of 

at least Aft, A <pj, and the LU-factorization of matrix Mi or the matrix inverse of matrix 
Mi; 

a fifth set of instructions embodied on the computer-readable medium which 
directs the data processor to generate, for an additional time moment j 5*1 , an LU- 
factorization of a matrix M, or a matrix inverse of matrix My, the matrix My being a 

function of at least A and Hj y ; and 

a sixth set of instructions embodied on the computer-readable medium which 
directs the data processor to generate, for an additional time moment j 5*1, a vector Ny as a 
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function of at least AjJ, A g)j , and the LU- factorization or matrix M, or the matrix inverse 
of matrix M j9 the vector N 7 having estimates of the floating ambiguities. 

24. The computer program product of Claim 23 wherein the third set of 
instructions directs the data processor to generate an LU-factorization for a matrix 

comprising a form equivalent to (g T Pig), where: 

• the matrix G has In rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprises annxn zero matrix and the other sub-matrix 
comprising an«x« identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

• the matrix Pi has 2n rows, 2n columns, and a form which comprises a matrix 

equivalent to ^Rf 1 - Rf 1 Qi(q^ Rf 1 Q x +?S 1 )~ 1 Qj r Rf 1 j where the matrix Ri 

is a weighting matrix, where the matrix Ri" 1 comprises an inverse of matrix R i5 where 
the matrix Qi has an upper sub-matrix and a lower sub-matrix, one of the sub- 
matrices of Qi comprising matrix /// 7 and the other of the sub-matrices of Qi 

comprising the matrix product A Hj y 9 and wherein the matrix Qi T comprises the 
transpose of matrix Qi, and where the quantity qSk is a zero matrix when the distance 
between the first and second navigation receivers is unconstrained and where q may 
be a non-zero weighting parameter and may be a non-zero matrix when the 
distance between the first and second navigation receivers is constrained. 

25. The computer program product of Claim 24 wherein the fourth set of 
instructions directs the data processor to generate vector Ni to comprise a vector having a 

form equivalent to M f 1 (g T Pj pi + ggi ], where the matrix Mr 1 comprises an inverse 
of matrix of matrix Mi, and where the vector fi\ comprises the vector [A/\,A<pi ] T , and 
where the quantity qgk is a zero vector when the distance between the first and second 
navigation receivers is unconstrained and where q may be non-zero and may be a non- 
zero vector when the distance between the first and second navigation receivers is 
constrained. 
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26. The computer program product of Claim 23 wherein the fifth set of 
instructions directs the data processor to generate an LU-factorization for a matrix 

T 

comprising a form equivalent to My = M^-.j + G PyG , where: 

• Mj.j comprises the matrix Mi of step (c) when j = 2 and comprises the matrix My of 
step (e) for the j-l time moment when j > 2, 

• the matrix G has 2n rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprising annxn zero matrix and the other sub-matrix 
comprising znnxn identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

• the matrix Py has 2n rows, 2n columns, and a form which comprise a matrix 



r 

equivalent to 



Ry 1 - R/ Qy(Qj R} 1 Qj +qSj)~ 1 Q r j R}^ where the matrix 

Ry is a weighting matrix, where the matrix Ry -1 comprises an inverse of matrix Ry, 
where the matrix Qy has an upper sub-matrix and a lower sub-matrix, one of the sub- 
matrices of Qy comprising matrix Hj* and the other of the sub-matrices of Qy 

comprising the matrix product A Hj*, and wherein the matrix Qy T comprises the 
transpose of matrix Qy, and where the quantity qSy is a zero matrix when the distance 
between the first and second navigation receivers is unconstrained and where q may 
be a non-zero weighting parameter and Sy may be a non-zero matrix when the distance 
between the first and second navigation receivers is constrained. 

27. The computer program product of Claim 26 wherein the sixth set of 
instructions directs the data processor to generate a vector N, to comprise a vector having 

a form equivalent to ^j-i + M^ 1 |G T Py ( fly - GNy^ )+ qgj\, where the matrix 

My" 1 comprises an inverse of matrix of matrix My, where the vector p,j comprises the 
vector [Ayj,dq>j] T 9 and where the vector Ny_y comprises the vector Ni generated by step 
(d) when j = 2 and comprises the vector Ny_y generated by step (f) for the j-l time moment 
when j > 2, and where the quantity qgj is a zero vector when the distance between the 
first and second navigation receivers is unconstrained and where q may be non-zero and 
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gj may be a non-zero vector when the distance between the first and second navigation 
10 receivers is constrained. 



10 



15 



10 



28. The computer program product of Claim 27 wherein the distance between the 
first and second navigation receivers is constrained to a distance Lrb, and wherein the 
fifth set of instructions directs the data processor to generate matrix Sj in a form 
equivalent to: 



1- 



JWJ 



0 3 xl 

0 



where r y is a vector comprising estimates of the three coordinates of the baseline vector 
for the j-th time moment and a zero as fourth vector component, where r/ is the vector 
transpose of r,, where I 3 is the 3 -by- 3 identity matrix, where Oi X 3 is a row vector of three 
zeros, and where C>3 X i is a column vector of three zeros; and 

wherein the sixth set of instructions directs the data processor to generate vector g, 
for the j-th time moment in a form equivalent to: 



gj =G T R7 1 Q,-(Q/RJ 1 Q,. + qSjy 1 hj, 



where: 



J\\J 



29. The computer program product of Claim 28 wherein the fifth set of 
instructions comprises: 

a seventh set of instructions that direct the data processor to generate an LU- 
factorization of matrix Mj.j in a form equivalent to L/_i L>j-\ T wherein L/_i is a low- 
triangular matrix and L/_i T is the transpose of L/_i ; 

an eighth set of instructions that direct the data processor to generate a 
factorization of G T P, G in a form equivalent to T, T/ = G T P, G, where T/ is the 
transpose of T, ; and 

a ninth set of instructions that direct the data processor to generate an LU- 
factorization of matrix M, in a form equivalent to L y luf from a plurality ft of rank-one 



75 



modifications of matrix L/_i, each rank-one modification being based on a respective 
column of matrix T y , where n is the number of rows in matrix M y . 

30. The computer program product of Claim 29 wherein the eighth set of 
instructions directs the data processor to generate matrix T y from a Cholesky factorization 
ofG T P,G. 



31 . The computer program product of Claim 29 wherein weighting matrix R 7 has 
a form equivalent to: 



R j = 



R] o 
0 Rj J 



, where R Y and R^ are weighting matrices; 



wherein R Y and R^ are related to a common weighting matrix W and scaling 

parameters a y and cr^ as follows (r 7 )" 1 =~y w > and (r^)~ = — y W ; and 

(7 y 

wherein the eighth set of instructions comprises: 

instructions that direct the data processor to generate a scalar b in a form 

CJyAQps 2 

equivalent to: b = — 2 2 2 9 w ^ ere ^GPS ls the wavelength of the satellite signals, 

CTy +^GPS cr <p 

instructions that direct the data processor to generate a matrix H in a form 
equivalent to H = W 1/2 HjJ , 

instructions that direct the data processor to generate a Householder matrix Shh 
for matrix H , and 

instructions that direct the data processor to generate matrix T) in a form 
equivalent to: 



T/= — W 1/2 Shh 



V A GPS J 



0 4x ( n _ 4 ) 



>(„-4)x4 I I( w -4)x(«-4) 
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32. The method of Claim 29 wherein weighting matrix R, is applied to a case 
where there is a first group of satellite signals having carrier signals in a first wavelength 
band and a second group of satellite signals having a carrier signals in a second 
wavelength band, the weighting frequency having a form equivalent to: 



r, j 



o 



0 Rj 



, where R 7 and R^ are weighting matrices; 



earner 



wherein R Y and R^ are related to a common weighting matrix W, the carrier 
wavelengths of the first group of signals as represented by matrix , the i 
wavelengths of the second group of signals as represented by matrix 

A <2) 

, the center 

wavelength of the first band as represented by , the center wavelength of the first band 
as represented by X 2 , and scaling parameters <j y and , as follows: 
1 



w 



o 



o 



nxn 



nxn 



w 



kr- 



o 



w 



(1) 



o 



nxn 



nxn 



J_ W <2) 



'<P 



15 



where 



W 0) = 1 A 0)^ A (l) w m = J_ A (2 >^A (2) 



wherein step (h) of generating matrix T 7 - comprises the steps of: 
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generating a scalar b in a form equivalent to: b — 



2 i2 5 



where >lj is the wavelength of a first group of satellite signals and is the wavelength 
of a second group of satellite signals, 

generating a matrix H in a form equivalent to H = W 1/2 H^ , 

generating a Householder matrix Shh for matrix H , and 

All I 



generating matrix 1} in a form equivalent to: Tj = 



1 



•<p 



A21 



'nxn 



A22 



where sub-matrixces All, A21, and A22 are as follows: 

yjl-b/Xi I 4 | 0 4x ( n _ 4 ) 



All 



= (w«) 



1 
2 S 



HH 



| 

°(n-4)x4 I '(11-4) 



A21 



= (w< 2 >) 2 



2 S 



HH 



^X^l-b/rf 



*4 I °4x(n-4) 



I 



°(n-4)x4 I °(n-4)x(n-4) 



, and 



A22 



= (w< 2 )) 2 : 



'HH 



\-b/X*-b/X* 

" — 2 14 ' °4x(n-4) 



\-b/, 



o (n _ 



(n-4)x4 



I(n-4) 



33. A computer program product for directing a data processor to estimate a set of 
floating ambiguities associated with a set of phase measurements of a plurality n of 
satellite carrier signals made by a first navigation receiver (B) and a second navigation 
receiver (R) separated by a distance, wherein a baseline vector (x°,y°,z 0 ) relates the 
position of the second receiver to the first receiver, each satellite carrier signal being 
transmitted by a satellite and having a wavelength, wherein each receiver has a time clock 
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for referencing its measurements and wherein any difference between the time clocks may 
be represented by an offset, the process receiving, for a plurality of two or more time 
moments j 9 the following inputs: 

a vector yf representative of a plurality of pseudo-ranges measured by the first 
navigation receiver (B) and corresponding to the plurality of satellite carrier signals, 

a vector yf representative of a plurality of pseudo-ranges measured by the second 
navigation receiver (R) and corresponding to the plurality of satellite carrier signals, 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the first navigation receiver (B), 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the second navigation receiver (R), 

a vector <Pj representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the first navigation receiver (B), 

a vector <pj representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the second navigation receiver (R), 

a geometric Jacobian matrix Hf whose matrix elements are representative of the 
changes in the distances between the satellites and one of the receivers that would be 
caused by changes in that receiver's position and time clock offset, the computer program 
product comprising: 

a computer-readable medium; 

a first set of instructions embodied on the computer-readable medium which 
directs the data processor to generate, for each time moment y, a vector A# of a plurality 
of range residuals of pseudo-range measurements made by the first and second navigation 
receivers in the form of: Aft = (yf -yf)- (Df - Df ); 

a second set of instructions embodied on the computer-readable medium which 
directs the data processor to generate, for each time moment j 9 a vector A(pj of a plurality 
of phase residuals of full phase measurements made by the first and second navigation 
receivers in the form of: A q>j = (<pf - <pf) - A 1 (Df - Df ), where A 1 is a diagonal 
matrix comprising the inverse wavelengths of the satellites; 
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a third set of instructions embodied on the computer-readable medium which 
directs the data processor to generate an LU- factorization of a matrix M or a matrix 

inverse of matrix M, the matrix M being a function of at least A and H/?, for index k 
of //^covering at least two of the time moments j ; and 

a fourth set of instructions embodied on the computer-readable medium which 
directs the data processor to generate a vector N of estimated floating ambiguities as a 

function of at least the set of range residuals A% , the set of phase residuals A (fa, and the 
LU- factorization of matrix M or the matrix inverse of matrix M. 

34. The computer program product of Claim 33 wherein the third set of 
instructions directs the data processor to generate matrix M in a form equivalent to the 
summation 

£(G T P*G) , 

_*=1 

where: 

• the matrix G has 2n rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprises annxn zero matrix and the other sub-matrix 
comprising annxw identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

• the matrix V k has 2n rows, 2n columns, and a form which comprises a matrix 

equivalent to = R* 1 - R* 1 qJqJ R* 1 Q k + qS k Y Qj R* 1 , where the 
matrix R* is a weighting matrix, where the matrix R*" 1 comprises an inverse of matrix 
Rjt, where the matrix Qk has an upper sub-matrix and a lower sub-matrix, one of the 
sub-matrices of comprising matrix and the other of the sub-matrices of 

comprising the matrix product A~ and wherein the matrix Q, T comprises the 
transpose of matrix Q*, and where the quantity qSk is a zero matrix when the distance 
between the first and second navigation receivers is unconstrained and where q may 
be a non-zero weighting parameter and S* may be a non-zero matrix when the distance 
between the first and second navigation receivers is constrained. 
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35. The computer program product of Claim 34 wherein the fourth set of 
instructions directs the data processor to generate matrix N in a form equivalent to: 



£(G T P*ji*+<7g k ) 

k=\ 



N = M _1 x 

where the matrix M" 1 comprises an inverse of matrix of matrix M, where the vector fa 
comprises the vector [Ay k ,A<p k ] T , and where the quantity qgk is a zero vector when the 
distance between the first and second navigation receivers is unconstrained and where q 
may be non-zero and gk may be a non-zero vector when the distance between the first and 
second navigation receivers is constrained. 



36. The computer program product of Claim 35 wherein the distance between the 
first and second navigation receivers is constrained to a distance Lrb, wherein the third set 
of instructions directs the data processor to generate matrix Sk for the k-th time moment 
in a form equivalent to: 



1- 



i 3 o 3xl 

IMA°i*3 0 . 



+ ±M r r k r k T 



Ik I 



where Y k is a vector comprising estimates of the three coordinates of the baseline vector at 
the k-th time moment, and a zero as fourth component, where r* T is the vector transpose 
of where I 3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three zeros, and 
where C>3 X i is a column vector of three zeros; and 

wherein step (d) generates vector gk for the k-th time moment in a form equivalent 

to: 



g k =G T R jt 1 Q t (Q jfc T R jt 1 Q A + 9 S ifc r 1 h i t, 



where: 



1- 



37. An apparatus for estimating a set of floating ambiguities associated with a set 
of phase measurements of a plurality n of satellite carrier signals made by a first 
navigation receiver (B) and a second navigation receiver (R) separated by a distance, 
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wherein a baseline vector (x°,y°,z°) relates the position of the second receiver to the first 
5 receiver, each satellite carrier signal being transmitted by a satellite and having a 

wavelength, wherein each receiver has a time clock for referencing its measurements and 
wherein any difference between the time clocks may be represented by an offset, said 
apparatus receiving, for a plurality of two or more time moments j, the following inputs: 

a vector yf representative of a plurality of pseudo-ranges measured by the first 
10 navigation receiver (B) and corresponding to the plurality of satellite carrier signals, 

a vector yf representative of a plurality of pseudo-ranges measured by the second 
navigation receiver (R) and corresponding to the plurality of satellite carrier signals, 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the first navigation receiver (B), 
15 a vector Dj representative of a plurality of estimated distances between the 

satellites and the second navigation receiver (R), 

a vector <Pj representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the first navigation receiver (B), 

a vector g>j representative of a plurality of full phase measurements of the 
20 satellite carrier signals measured by the second navigation receiver (R), 

a geometric Jacobian matrix Hj y whose matrix elements are representative of the 
changes in the distances between the satellites and one of the receivers that would be 
caused by changes in that receiver's position and time clock offset, 

said apparatus comprising: 

25 (a) means for generating, for each time moment j, a vector A# of a plurality of 

range residuals of pseudo-range measurements made by the first and second navigation 
receivers in the form of: A# = {yf - yf) - {Df - Df ); 

(b) means for generating, for each time moment j , a vector A q>j of a plurality of 
phase residuals of full phase measurements made by the first and second navigation 

R B — 1 — 1 

30 receivers in the form of: A q)j = (tpj - q)j ) — A (Dj - Dj ), where A is a diagonal 
matrix comprising the inverse wavelengths of the satellites; 
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(c) means for generating, for time moment j= 1, an LU-factorization of a matrix 

. -l 

Mi or a matrix inverse of matrix Mi, the matrix Mi being a function of at least A and 

35 (d) means for generating, for time moment j = 1, a vector Ni as a function of at 

least Aft, A and the LU-factorization of matrix Mi or the matrix inverse of matrix 



M 



i; 



(e) means for generating, for an additional time moment j 5*1, an LU-factorization 
of a matrix My or a matrix inverse of matrix My, the matrix My being a function of at least 

40 A~* and H/; and 

(f) means for generating, for an additional time moment j 5*1, a vector Ny as a 
function of at least Afi 9 A q>j , and the LU-factorization or matrix My or the matrix inverse 
of matrix My, the vector Ny having estimates of the floating ambiguities. 

38. The apparatus of Claim 37 wherein means (c) comprises means for generating 
an LU-factorization for a matrix comprising a form equivalent to (g T Pig), where: 

• the matrix G has 2n rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprises annxn zero matrix and the other sub-matrix 

5 comprising annxn identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

• the matrix Pi has 2n rows, In columns, and a form which comprises a matrix 

equivalent to Rf 1 - Rf 1 Q { (q^ Rf 1 Q 2 + qSi ^Qi Rf 1 j where the matrix Ri 

is a weighting matrix, where the matrix Ri" 1 comprises an inverse of matrix R u where 
10 the matrix Qi has an upper sub-matrix and a lower sub-matrix, one of the sub- 

matrices of Qi comprising matrix //j 7 and the other of the sub-matrices of Qi 

comprising the matrix product A Hj\ and wherein the matrix Qi T comprises the 
transpose of matrix Qi, and where the quantity qSk is a zero matrix when the distance 
between the first and second navigation receivers is unconstrained and where q may 
15 be a non-zero weighting parameter and S# may be a non-zero matrix when the 

distance between the first and second navigation receivers is constrained. 
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39. The apparatus of Claim 38 wherein means (d) comprises means for generating 

vector Ni to comprise a vector having a form equivalent to M j" 1 (g T Pj pi + qgi Y where 
the matrix Mi" 1 comprises an inverse of matrix of matrix Mi, and where the vector fi\ 
comprises the vector [A^i, Aq>\ ] T , and where the quantity qgk is a zero vector when the 
distance between the first and second navigation receivers is unconstrained and where q 
may be non-zero and gk may be a non-zero vector when the distance between the first and 
second navigation receivers is constrained. 

40. The apparatus of Claim 37 wherein means (e) comprises means for generating 
an LU- factorization for a matrix comprising a form equivalent to 

M y = Mj_ x + G T P y -G , where: 

• M/_/ comprises the matrix Mi of generated by means (c) when j = 2 and comprises the 
matrix My generated by means (e) for the y-1 time moment when j > 2, 

• the matrix G has In rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprising annxn zero matrix and the other sub-matrix 
comprising anwxn identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

• the matrix P, has 2n rows, 2n columns, and a form which comprise a matrix 

equivalent to ^RJ 1 - Ry 1 Qy(<2y Ry 1 Qy + ?Sy)~ 1 Qy RJ*j where the matrix 

Kj is a weighting matrix, where the matrix R/ 1 comprises an inverse of matrix R y , 
where the matrix Q/ has an upper sub-matrix and a lower sub-matrix, one of the sub- 
matrices of Qj comprising matrix H/and the other of the sub-matrices of Q, 

comprising the matrix product A Hj*, and wherein the matrix Q/ comprises the 
transpose of matrix Q/, and where the quantity qSy is a zero matrix when the distance 
between the first and second navigation receivers is unconstrained and where q may 
be a non-zero weighting parameter and S, may be a non-zero matrix when the distance 
between the first and second navigation receivers is constrained. 
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41. The apparatus of Claim 40 wherein means (f) comprises means for generating 
vector Ny to comprise a vector having a form equivalent to 

Ny_! + My 1 [G T Py ( fij - GNj_ x )+ qgj J, where the matrix M/ 1 comprises an 

inverse of matrix of matrix My, where the vector fij comprises the vector [Ay),J(pj ] T , and 
where the vector N^y comprises the vector Ni generated by means (d) when j = 2 and 
comprises the vector N y _/ generated by means (f) for the j-l time moment when j > 2, and 
where the quantity qgj is a zero vector when the distance between the first and second 
navigation receivers is unconstrained and where q may be non-zero and gj may be a non- 
zero vector when the distance between the first and second navigation receivers is 
constrained. 



42. The apparatus of claim 39 wherein the distance between the first and second 
navigation receivers is constrained to a distance Lrb, wherein means (c) comprises means 
for generating matrix Si in a form equivalent to: 



1- 



ll r l|| A°lx3 



°3xl 
0 



, L RB r T 

+ mr r i ri 
Fill 



where l*i is a vector comprising estimates of the three coordinates of the baseline vector 
for the time moment j=l and a zero as fourth component, where l"i T is the vector 
transpose of where I 3 is the 3-by-3 identity matrix, where Oi x3 is a row vector of three 
zeros, and where 0 3x i is a column vector of three zeros; and 

wherein means (d) comprises means for generating vector gi for the time moment 
j=l in a form equivalent to: 

gl ^gTr^q^q^r^Q! +< 7 s 1 r 1 h 1 , 

where: 



hi = 



1- 



'RB 



r l 



43. The apparatus of Claim 41 wherein the distance between the first and second 
navigation receivers is constrained to a distance Lrb, wherein means (e) comprises means 
for generating matrix Sj in a form equivalent to: 
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-RB 



h 
\P\x3 



°3xl 
0 



where Yj is a vector comprising estimates of the three coordinates of the baseline vector 
for the j-th time moment and a zero as fourth vector component, where r/ is the vector 
transpose of Ty, where I3 is the 3-by-3 identity matrix, where Oi x3 is a row vector of three 
zeros, and where 03 X i is a column vector of three zeros; and 

wherein means (f) comprises means for generating vector g, for the j-th time 
moment in a form equivalent to: 

gj ^G^fQjiQ/RfQj +< ? S y )- 1 h / , 

where: 



1- 



'RB 



J\\J 



44. The apparatus of Claim 40 wherein means (e) for generating the LU- 
factorization of My comprises: 

(g) means for generating an LU- factorization of matrix Mj.j in a form equivalent 
to Ly_i L/_i r wherein L y _i is a low-triangular matrix and L/_i r is the 
transpose of L/_i; 

(h) means for generating a factorization of G T Py G in a form equivalent to Ty Ty T = 
G T Py G, where Ty T is the transpose of Ty ; and 

(i) means for generating an LU-factorization of matrix My in a form equivalent to 
Ly Ly r from a plurality n of rank-one modifications of matrix Ly_j, each rank-one 
modification being based on a respective column of matrix Ty, where n is the number of 
rows in matrix My. 

45. The apparatus of Claim 44 wherein means (h) generates matrix Ty from a 
Cholesky factorization of G T Py G. 
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46. The apparatus of Claim 44 wherein weighting matrix R, has a form equivalent 



to: 



R j = 



RZ 0 
0 R? 



, where R Y and are weighting matrices; 



wherein R Y and R^ are related to a common weighting matrix W and scaling 
parameters a y and as follows (r y }^ =^tW, and (r^)" 1 = — y W, 



wherein means (h) of generating matrix T, comprises: 



& Y GPS 

means for generating a scalar b in a form equivalent to: b = — — , where 



2 ^ il__„2 
r 



<T V +%GPS (y <p 



&GPS ls ^ wavelength of the satellite signals, 



means for generating a matrix H in a form equivalent to H = W 1/2 H^ , 

means for generating a Householder matrix S H h for matrix H , and 
means for generating matrix T) in a form equivalent to: 



T> = — W 1/2 Shh 



1- 



0(„_ 4)x4 



I4 I o 4x ( W _ 4) 
I 

I l (n-4)x(n-4) 



47. The method of Claim 44 wherein weighting matrix R, is applied to a case 
where there is a first group of satellite signals having carrier signals in a first wavelength 
band and a second group of satellite signals having a carrier signals in a second 
wavelength band, the weighting frequency having a form equivalent to: 



R j = 



R) 0 
0 Rj 



, where R Y and R^ are weighting matrices; 



wherein R Y and R^ are related to a common weighting matrix W, the carrier 
wavelengths of the first group of signals as represented by matrix .the 
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wavelengths of the second group of signals as represented by matrix , the center 
wavelength of the first band as represented by , the center wavelength of the first band 

as represented by X 2 , and scaling parameters cj y and cr^ , as follows: 

1 



for- 



W I o 



for- 



o 

1 

2 

o 



nxn 



W 



(1) 



nxn 



nxn 



w 



o 



nxn 



J_ W (2, 



<P 



15 



20 



where W^=\ A« WA^ W™ = \ A (2) PTA (2) 
wherein step (h) of generating matrix Tj comprises the steps of: 



generating a scalar b in a form equivalent to: ^ 



where /Lj is the wavelength of a first group of satellite signals and A 2 is the wavelength 
of a second group of satellite signals, 

generating a matrix H in a form equivalent to H = W 1/2 H^ , 

generating a Householder matrix S H h for matrix H , and 

All I O t 



1 



generating matrix T, in a form equivalent to: Tj = 



'nxn 



A21 | A22 



where sub-matrixces All, A21, and A22 are as follows: 
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All 



= (W< 1 ))2 



S HH 



jl-b/tf I 4 I 0 4x(n _ 4) 
I 

°(n-4)x4 I J (n-4) 



A21 



= (\V( 2 >)2 



2 S 



HH 



°(n-4)x4 



I4 I °4x(n-4) 



°(n-4)x(n-4) 



, and 



A22 



= (w< 2 >)*: 



•HH 



\\-b/X* -VX* 



I °4x(n-4) 
°(n-4)x4 I I(n-4) 



48. An apparatus for estimating a set of floating ambiguities associated with a set 
of phase measurements of a plurality n of satellite carrier signals made by a first 
navigation receiver (B) and a second navigation receiver (R) separated by a distance, 
wherein a baseline vector (x° f y 0 ,z°) relates the position of the second receiver to the first 
receiver, each satellite carrier signal being transmitted by a satellite and having a 
wavelength, wherein each receiver has a time clock for referencing its measurements and 
wherein any difference between the time clocks may be represented by an offset, said 
apparatus receiving, for a plurality of two or more time moments j 9 the following inputs: 

a vector yf representative of a plurality of pseudo-ranges measured by the first 
navigation receiver (B) and corresponding to the plurality of satellite carrier signals, 

a vector yf representative of a plurality of pseudo-ranges measured by the second 
navigation receiver (R) and corresponding to the plurality of satellite carrier signals, 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the first navigation receiver (B), 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the second navigation receiver (R), 



89 



a vector q)j representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the first navigation receiver (B), 

a vector q)j representative of a plurality of full phase measurements of the 
20 satellite carrier signals measured by the second navigation receiver (R), 

a geometric Jacobian matrix Hj y whose matrix elements are representative of the 
changes in the distances between the satellites and one of the receivers that would be 
caused by changes in that receiver's position and time clock offset, 

said apparatus comprising: 

25 (a) means for generating, for each time moment j 9 a vector Aft of a plurality of 

range residuals of pseudo-range measurements made by the first and second navigation 
receivers in the form of: Aft = (yf - yf) - (Df - D* ), said means generating a set of 
range residuals A% , fc=l, „,y; 

(b) means for generating, for each time moment j, a vector A q>j of a plurality of 
30 phase residuals of fxill phase measurements made by the first and second navigation 

receivers in the form of: A = (ty - <Pj ) — A (Dj -Dj ), where A is a diagonal 
matrix comprising the inverse wavelengths of the satellites, said means generating a set of 
phase residuals Aq>k , fc=l, ..J\ 

(c) means for generating an LU-factorization of a matrix M or a matrix inverse of 
35 matrix M, the matrix M being a function of at least A and for index k of Hk 

covering at least two of the time moments j ; 

(d) means for generating a vector N of estimated floating ambiguities as a function 
of at least the set of range residuals Ay k , the set of phase residuals A qk, and the LU- 
factorization of matrix M or the matrix inverse of matrix M. 



49. The apparatus of Claim 48 wherein means (c) comprises means for generating 
matrix M in a form equivalent to the summation 
where: 



Z(G T P*G) 
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• the matrix G has 2n rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprises an n x n zero matrix and the other sub-matrix 
comprising an/ix« identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

• the matrix P* has 2n rows, In columns, and a form which comprises a matrix 

equivalent to = R* 1 - R* 1 Q k (q J R* 1 Q k + qS k Qj R* 1 , where the 
matrix R^ is a weighting matrix, where the matrix R* -1 comprises an inverse of matrix 
R*, where the matrix has an upper sub-matrix and a lower sub-matrix, one of the 
sub-matrices of Q k comprising matrix and the other of the sub-matrices of 
comprising the matrix product A H k \ and wherein the matrix Q* T comprises the 
transpose of matrix Q*, and where the quantity qSk is a zero matrix when the distance 
between the first and second navigation receivers is unconstrained and where q may 
be a non-zero weighting parameter and S* may be a non-zero matrix when the distance 
between the first and second navigation receivers is constrained. 

50. The apparatus of Claim 49 wherein means (d) comprises means for generating 
matrix N in a form equivalent to: 

X(G T P*^+<7g k ) , 

k=l 

where the matrix M" 1 comprises an inverse of matrix of matrix M, where the vector fi k 
comprises the vector [Ay k ,/i<p k ] T , and where the quantity qgk is a zero vector when the 
distance between the first and second navigation receivers is unconstrained and where q 
may be non-zero and may be a non-zero vector when the distance between the first and 
second navigation receivers is constrained. 

5 1 . The apparatus of claim 50 wherein the distance between the first and second 
navigation receivers is constrained to a distance Lrb, wherein means (c) comprises means 
for generating matrix S k for the k-th time moment in a form equivalent to: 

c, _fi L RB Y 13 °3xQ Lrb_ T 

*t raJl<>.*3 0 j + w r * r * 



N = M _1 x 
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where r k is a vector comprising estimates of the three coordinates of the baseline vector at 
the k-th time moment, and a zero as fourth component, where r* T is the vector transpose 
of r*, where I3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three zeros, and 
where 03 X i is a column vector of three zeros; and 

wherein means (d) generates vector g k for the k-th time moment in a form 
equivalent to: 

g k = G T R* 1 Q*(Q* T R* 1 Q* + qS k y 1 h k , 



where: 



52. (Claims for Japan) A computer program to be installed in a computer for 
controlling the computer to perform the process of estimating a set of floating ambiguities 
associated with a set of phase measurements of a plurality n of satellite carrier signals 
made by a first navigation receiver (B) and a second navigation receiver (R) separated by 
a distance, wherein a baseline vector (x°,y°,z°) relates the position of the second receiver 
to the first receiver, each satellite carrier signal being transmitted by a satellite and having 
a wavelength, wherein each receiver has a time clock for referencing its measurements 
and wherein any difference between the time clocks may be represented by an offset, said 
method receiving, for a plurality of two or more time moments j\ the following inputs: 

a vector yf representative of a plurality of pseudo-ranges measured by the first 
navigation receiver (B) and corresponding to the plurality of satellite carrier signals, 

a vector representative of a plurality of pseudo-ranges measured by the second 
navigation receiver (R) and corresponding to the plurality of satellite carrier signals, 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the first navigation receiver (B), 

a vector Dj representative of a plurality of estimated distances between the 



satellites and the second navigation receiver (R), 

a vector q>j representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the first navigation receiver (B), 
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20 a vector <Pj representative of a plurality of full phase measurements of the 

satellite carrier signals measured by the second navigation receiver (R), 

a geometric Jacobian matrix Hj y whose matrix elements are representative of the 
changes in the distances between the satellites and one of the receivers that would be 
caused by changes in that receiver's position and time clock offset, said process 

25 comprising: 

(a) generating, for each time moment J, a vector A^- of a plurality of range 
residuals of pseudo-range measurements made by the first and second navigation 
receivers in the form of: A% = ~yf)~ (Dj* - D* ); 

(b) generating, for each time moment j 9 a vector A <pj of a plurality of phase 

30 residuals of full phase measurements made by the first and second navigation receivers in 
the form of: 

A_<pj = (^ R - <p?) - A"' iDf - Df ), 
where A is a diagonal matrix comprising the inverse wavelengths of the satellites; 

(c) generating, for time moment j = 1, an LU- factorization of a matrix Mi or a 
35 matrix inverse of matrix Mi, the matrix Mi being a function of at least A and H y ; 

(d) generating, for time moment j = 1, a vector Ni as a function of at least A^, 

A q)j, and the LU-factorization of matrix Mi or the matrix inverse of matrix Mi; 

(e) generating, for an additional time moment j 5*1, an LU-factorization of a 

. -l 

matrix M, or a matrix inverse of matrix M 7 , the matrix My being a function of at least A , 
40 Hj y and an instance of matrix M generated for a different time moment; and 

(f) generating, for an additional time moment j 5*1 , a vector Ny as a function of at 
least A#, A <pj , and the LU-factorization or matrix My or the matrix inverse of matrix My, 
the vector Ny having estimates of the floating ambiguities. 

53. The computer program of Claim 52 wherein step (c) of the process comprises 

generating an LU-factorization for a matrix comprising a form equivalent to (g^jg), 
where: 
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the matrix G has In rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprises annxn zero matrix and the other sub-matrix 
comprising an n x n identity matrix, 

the matrix G T comprises the transpose matrix of matrix G, and 

the matrix Pi has 2n rows, 2n columns, and a form which comprises a matrix 

equivalent to ^Rf 1 - Rf 1 Q^Q^ Rf 1 Qi + ?S 1 )~ 1 Qj r Rf !> where the matrix Ri 

is a weighting matrix, where the matrix Ri -1 comprises an inverse of matrix Ri, where 
the matrix Qi has an upper sub-matrix and a lower sub-matrix, one of the sub- 
matrices of Q, comprising matrix JZ/and the other of the sub-matrices of Q, 

comprising the matrix product A Hj\ and wherein the matrix Qi T comprises the 
transpose of matrix Qi, and where the quantity qSk is a zero matrix when the distance 
between the first and second navigation receivers is unconstrained and where q may 
be a non-zero weighting parameter and Sjt may be a non-zero matrix when the 
distance between the first and second navigation receivers is constrained. 



54. The computer program of Claim 53 wherein step (d) of the process comprises the 
step of generating vector Ni to comprise a vector having a form equivalent to 

Mj- 1 (G T P 1 

Pi + )> where the matrix M, 1 comprises an inverse of matrix of matrix 
Mi, and where the vector /ii comprises the vector [Ay u A<p\ ] T , and where the quantity 
qgk is a zero vector when the distance between the first and second navigation receivers 
is unconstrained and where q may be non-zero and gk may be a non-zero vector when the 
distance between the first and second navigation receivers is constrained. 

55. The computer program of Claim 52 wherein step (e) of the process comprises 
generating an LU-factorization for a matrix comprising a form equivalent to 

My - My_! + G T P y G , where: 

• Mj.j comprises the matrix Mi of step (c) when j = 2 and comprises the matrix My of 
step (e) for the j-l time moment when j > 2, 
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• the matrix G has In rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
one of the sub-matrices comprising annxn zero matrix and the other sub-matrix 
comprising annxn identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

10 • the matrix P, has 2n rows, 2n columns, and a form which comprise a matrix 

equivalent to ^RJ 1 - R} 1 Q,(qJ R} 1 Q y + #S y )~ 1 Qj where the matrix 

R, is a weighting matrix, where the matrix R, -1 comprises an inverse of matrix R,, 
where the matrix Qy has an upper sub-matrix and a lower sub-matrix, one of the sub- 
matrices of Q/ comprising matrix Hj* and the other of the sub-matrices of Qy 

1 5 comprising the matrix product A Hj\ and wherein the matrix Q/ T comprises the 

transpose of matrix Q,, and where the quantity qSy is a zero matrix when the distance 
between the first and second navigation receivers is unconstrained and where q may 
be a non-zero weighting parameter and Sy may be a non-zero matrix when the distance 
between the first and second navigation receivers is constrained. 

20 

56. The computer program of Claim 55 wherein step (f) of the process comprises the 
step of generating vector Ny to comprise a vector having a form equivalent to 

Ny_i 4- My 1 [g T P/ ( \ij — GNy_i )+ qgj j, where the matrix My" 1 comprises an 

inverse of matrix of matrix My, where the vector jiy comprises the vector [A#, 4q)j] T , and 
5 where the vector Ny_; comprises the vector Ni generated by step (d) when j = 2 and 

comprises the vector Ny_y generated by step (f) for the y-1 time moment when j > 2, and 
where the quantity qgj is a zero vector when the distance between the first and second 
navigation receivers is unconstrained and where q may be non-zero and gj may be a non- 
zero vector when the distance between the first and second navigation receivers is 
10 constrained. 

57. The computer program of claim 56 wherein the distance between the first and 
second navigation receivers is constrained to a distance Lrb, wherein step (c) of the 
process generates matrix Si in a form equivalent to: 
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I hlA°ix3 0 J ||ri|| 1 
where 1*1 is a vector comprising estimates of the three coordinates of the baseline vector 
for the time moment j=l and a zero as fourth component, where l*i T is the vector 
transpose of Ti, where I3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three 
zeros, and where Chxi is a column vector of three zeros; and 

wherein step (d) generates vector gi for the time moment j=l in a form equivalent 

to: 

gl = G T Rf 1 Q 1 (Q 1 T Rf 1 Q 1 +qrS 1 r 1 h l5 

where: 

( T \ 
h 1= r,. 

I Fill J 



58. The computer program of claim 56 wherein the distance between the first and 
second navigation receivers is constrained to a distance Lrb, wherein step (e) of the 
process generates matrix S, in a form equivalent to: 

f V 



S/ = 



1- 



'RB 



JWJ 



1 



o 



3x1 



Oi x3 0 



where r y is a vector comprising estimates of the three coordinates of the baseline vector 
for the j-th time moment and a zero as fourth vector component, where r/ is the vector 
transpose of r y , where I3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three 
zeros, and where C>3 X i is a column vector of three zeros; and 

wherein step (f) generates vector g, for the j-th time moment in a form equivalent 

to: 

gj =G T R7 1 Q y (Q/R7 1 Q, +flS / r 1 h /f 

where: 



1- 



Ml; 
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59. The computer program of Claim 55 wherein the generation of the LU- 
factorization in step (e) of the process comprises the steps of: 

(g) generating an LU- factorization of matrix My.; in a form equivalent to 
L,_i L 7 _i T wherein L/_j is a low-triangular matrix and L/_i r is the transpose of Ly_i; 
5 (h) generating a factorization of G T Py G in a form equivalent to T) T y T = G T Py G, 

where Ty T is the transpose of Ty ; 

(i) generating an LU-factorization of matrix My in a form equivalent to Ly L 7 r from 
a plurality # of rank-one modifications of matrix Ly_i, each rank-one modification being 
based on a respective column of matrix Ty, where n is the number of rows in matrix My . 

60. The computer program of Claim 59 wherein step (h) of the process generates 
matrix Ty from a Cholesky factorization of G T Py G. 

61. The computer program of Claim 55 wherein step (d) of the process comprises 
generating a vector Bi to comprise a vector having a form equivalent to G T Pi Mi +c Lgi> 
where the vector fi\ comprises the vector [Ay u J<p\ ] T , and where the quantity qgj is a 
zero vector when the distance between the first and second navigation receivers is 

5 unconstrained and where q may be non-zero and gj may be a non-zero vector when the 
distance between the first and second navigation receivers is constrained; and 

wherein step (f) of the process further comprises generating, for each time moment 
j ?*l 9 a. vector By to comprise a matrix having a form equivalent to By_i + G T Py fly +qg/ 5 
where the vector fij comprises the vector [A#, A <pj ] T , and where the vector B y _y is the 

10 vector Bi generated by step (d) when j = 2 and comprises the vector generated by step (f) 
for the for the j-1 time moment when j > 2, and where the quantity qg, is a zero vector 
when the distance between the first and second navigation receivers is unconstrained and 
where q may be non-zero and g, may be a non-zero vector when the distance between the 
first and second navigation receivers is constrained; and 

15 wherein step (f) of the process further comprises generating vector Ny to comprise 

a vector having a form equivalent to Ny = [My]" 1 By, where the matrix My" 1 comprises an 
inverse of matrix of matrix M,. 
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62. The computer program of claim 61 wherein the distance between the first and 
second navigation receivers is constrained to a distance Lrb, wherein step (c) of the 
process generates matrix Si in a form equivalent to: 



S, = 



1- 



r li 



v i 3 <w 



^°lx3 



0 



+ ^. riri 7- 



where 1*1 is a vector comprising estimates of the three coordinates of the baseline vector 
for the time moment j=l and a zero as fourth component, where i*i T is the vector 
transpose of T\, where I 3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three 
zeros, and where 0 3x i is a column vector of three zeros; and 

wherein step (d) generates vector gi for the time moment j=l in a form equivalent 

to: 

gl =G T Rf 1 Q 1 (Q 1 T Rf 1 Q 1 +< ? S 1 r 1 h 1 , 

where: 



h, = 



1- 



L RB 

hi; 



63. The computer program of claim 61 wherein the distance between the first and 
second navigation receivers is constrained to a distance Lrb, wherein step (e) of the 
process generates matrix Sy in a form equivalent to: 



1- 



-RB 



JWJ 



I 3 

V°lx3 



0 3 xl 

0 



where Tj is a vector comprising estimates of the three coordinates of the baseline vector 
for the j-th time moment and a zero as fourth vector component, where r/ is the vector 
transpose of r y , where I 3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three 
zeros, and where 0 3x i is a column vector of three zeros; and 

wherein step (f) of the process generates vector g, for the j-th time moment in a 
form equivalent to: 

gj=G T Rj 1 Q J (Q J T R- J 1 Qj + qS J y 1 h j , 

where: 
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64. (Claims for Japan) A computer program to be installed in a computer for 
controlling the computer to perform the process of estimating a set of floating ambiguities 
associated with a set of phase measurements of a plurality n of satellite carrier signals 
made by a first navigation receiver (B) and a second navigation receiver (R), wherein a 
baseline vector (x 0 ,y°,z°) relates the position of the second receiver to the first receiver, 
each satellite carrier signal being transmitted by a satellite and having a wavelength, 
wherein each receiver has a time clock for referencing its measurements and wherein any 
difference between the time clocks may be represented by an offset, said method 
receiving, for a plurality of two or more time moments j, the following inputs for each 
time moment j: 

a vector yf representative of a plurality of pseudo-ranges measured by the first 
navigation receiver (B) and corresponding to the plurality of satellite carrier signals, 

a vector yf representative of a plurality of pseudo-ranges measured by the second 
navigation receiver (R) and corresponding to the plurality of satellite carrier signals, 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the first navigation receiver (B), 

a vector Dj representative of a plurality of estimated distances between the 
satellites and the second navigation receiver (R), 

a vector q)j representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the first navigation receiver (B), 

a vector <pj representative of a plurality of full phase measurements of the 
satellite carrier signals measured by the second navigation receiver (R), 

a geometric Jacobian matrix //} T whose matrix elements are representative of the 
changes in the distances between the satellites and one of the receivers that would be 
caused by changes in that receiver's position and time clock offset, the process 
comprising: 
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(a) generating, for each time moment j 9 a vector Ayj of a plurality of range 
residuals of pseudo-range measurements made by the first and second navigation 
receivers in the form of: Ayj = (yf - rf) - (Dj R - Dj B ), said step generating a set of 

30 range residuals Ay k , £=1, „,y; 

(b) generating, for each time moment j, a vector A <pj of a plurality of phase 
residuals of full phase measurements made by the first and second navigation receivers in 
the form of: A = (^ R - q>j B ) — A 1 (Dj R - D* ), where A 1 is a diagonal matrix 
comprising the inverse wavelengths of the satellites, said step generating a set of phase 

35 residuals A q> k , h=l , ..J ; 

(c) generating an LU- factorization of a matrix M or a matrix inverse of matrix M, 
the matrix M being a function of at least A and H k , for index k oiH k covering at 
least two of the time moments j ; 

(d) generating a vector N of estimated floating ambiguities as a function of at least 
40 the set of range residuals A% , the set of phase residuals A 0Jfc, and the LU- factorization of 

matrix M or the matrix inverse of matrix M. 



65. The computer program of Claim 64 wherein step (c) of the process comprises 
generating matrix M in a form equivalent to the summation 



Z(G T P*G) 



_k=\ 
where: 

• the matrix G has 2n rows, n columns, an upper sub-matrix, and a lower sub-matrix, 
5 one of the sub-matrices comprises annxn zero matrix and the other sub-matrix 

comprising annxn identity matrix, 

• the matrix G T comprises the transpose matrix of matrix G, and 

• the matrix P* has 2n rows, 2n columns, and a form which comprises a matrix 

equivalent to = R^ 1 - R^ 1 Q^(qJ R^ 1 Q k +qS k \ 1 Qj R^ 1 , where the 
10 matrix R* is a weighting matrix, where the matrix R^ 1 comprises an inverse of matrix 

R*, where the matrix has an upper sub-matrix and a lower sub-matrix, one of the 
sub-matrices of comprising matrix H k y and the other of the sub-matrices of Q# 
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comprising the matrix product A H k y 9 and wherein the matrix Q^ T comprises the 
transpose of matrix Q*, and where the quantity qSk is a zero matrix when the distance 
1 5 between the first and second navigation receivers is unconstrained and where q may 

be a non-zero weighting parameter and S* may be a non-zero matrix when the distance 
between the first and second navigation receivers is constrained. 

66. The computer program of Claim 65 wherein step (d) of the process comprises 
generating matrix N in a form equivalent to: 

j 



N = M~*x 



k=\ 



where the matrix 1VT 1 comprises an inverse of matrix of matrix M, where the vector /i k 
comprises the vector [Ay k , A <p k ] T , and where the quantity qgk is a zero vector when the 
distance between the first and second navigation receivers is unconstrained and where q 
may be non-zero and gk may be a non-zero vector when the distance between the first and 
second navigation receivers is constrained. 

67. The computer program of claim 66 wherein the distance between the first and 
second navigation receivers is constrained to a distance Lrb, wherein step (c) of the 
process generates matrix S k for the k-th time moment in a form equivalent to: 



l_ L RB 



\ ir * ii j 



^> lx3 0 j 



, L RB „ r1 T 

+ li ii r * r * 

WkW 



5 where r k is a vector comprising estimates of the three coordinates of the baseline vector at 
the k-th time moment, and a zero as fourth component, where r* T is the vector transpose 
of r*, where I 3 is the 3-by-3 identity matrix, where Oi X 3 is a row vector of three zeros, and 
where 0 3x i is a column vector of three zeros; and 

wherein step (d) of the process generates vector gk for the k-th time moment in a 
10 form equivalent to: 



where: 
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h t = 



.-ft*. 



68. The computer program of Claim 65 wherein at least one of the weighting 
matrices R k comprises an identity matrix multiplied by a scalar quantity. 
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